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Abstract 

Atomically resolved imaging and spectroscopic characteristics of graphene grown by chemical vapor deposition 
(CVD) on copper are investigated by means of scanning tunneling microscopy and spectroscopy (STM/STS). For 
CVD-grown graphene remaining on the copper substrate, the monolayer carbon structures exhibit ripples and appear 
strongly strained, with different regions exhibiting different lattice structures and electronic density of states (DOS). 
In particular, ridges appear along the boundaries of different lattice structures, which exhibit excess charging effects. 
Additionally, the large and non-uniform strain induces pseudo-magnetic field up to ~ 50 Tesla, as manifested by the 
DOS peaks at quantized energies that correspond to pseudo-magnetic field-induced integer and fractional Landau lev- 
els. In contrast, for graphene transferred from copper to Si02 substrates after the CVD growth, the average strain on 
the whole diminishes, so do the corresponding charging effects and pseudo-magnetic fields except for sample areas 
near topological defects. These findings suggest feasible nano-scale "strain engineering" of the electronic states of 
graphene by proper design of the substrates and growth conditions. 

Keywords: Scanning tunneling microscopy, Scanning tunneling spectroscopies, Chemical vapor deposition, 
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1. Introduction 

The superior physical properties of graphene (]]] 
and its compatibility with two-dimensional lithographic 
processes have stimulated intense research of graphene- 
based electronics for "beyond Si-CMOS" technol- 
ogy Jl S 13 S S H . One of the major challenges for 
realizing the graphene-based beyond Si-CMOS tech- 
nology is the fabrication of high-quality large-area 
graphene sheets and the retention of superior electronic 
characteristics of graphene. The apparent degradation 
of carrier mobility and significant variations in the elec- 
tronic characteristics when graphene comes in contact 
with various dielectrics Jgl [loll suggest significant sus- 
ceptibility of the single-layer carbon atoms to the sur- 
rounding environment. On the other hand, the strong 
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susceptibility of graphene to external influences also 
provides opportunities for engineering unique proper- 
ties of graphene. For instance, it has been theoretically 
proposed that a designed strain aligned along three main 
crystallographic directions induces strong gauge fields 
that effectively act as a uniform magnetic field on the 
Dirac electrons 1 1 ill 1 211 . In particular, for a finite doping 
level, the quantizing field results in an insulating bulk 
and a pair of counter-circulating edge states, which is 
similar to the case of topological insulators lfl^[l4l[l5ll . 
Moreover, strained superlattices have been shown the- 
oretically to open significant energy gaps in the elec- 
tronic spectrum of graphene 111 ill , which could be much 
more effective than current approaches to the bandgap 
engineering of graphene by either controlling the width 
of graphene nano-ribbons I m l 1711 or applying magnetic 
fields to bi-laver graphene Hl8LTl9l [20h . Recent scanning 
tunneling microscopic and spectroscopic (STM/STS) 
studies of graphene nano-bubbles on Pt(lll) single 
crystal substrates II21II have verified the theoretical pre- 
diction of strain-induced pseudo-magnetic fields, and 
significant field values in excess of 300 Tesla have been 
observed. 
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In the context of synthesis of large area graphene, var- 
ious growth techniques have been developed to date, in- 
cluding ultra-high-vacuum annealing to desorption of Si 
from SiC single crystal surfaces l22i 12311 . deposition of 
graphene oxide films from a liquid suspension followed 
by chemical reduction 12412511 . and chemical vapor de- 
position (CVD) on transition metals such as Ru |2£ 



W, Ni S Hi, Co |30|], Pt ||3il |32j] and Cu 
While samples as large as centimeter squared have been 
achieved, and characterizations of macroscopic physical 
properties such as Raman scattering and mobility mea- 
surements have been reported J28, 2§l 33h, most large- 
area samples appear to exhibit thicknesses ranging from 
one to approximately twelve layers. Moreover, few in- 
vestigations with the exception of Refs. 02 1 1 13411 have 
been made on correlating the local electronic properties 
with the microscopic structural variations of these large- 
area films. 



In this work we report scanning tunneling micro- 
scopic and spectroscopic (STM/STS) investigations 
of the microscopic electronic and structural proper- 
ties of CVD-grown graphene on copper and also on 
SiC»2 after transferred from copper. We find that 
the CVD grown graphene is strongly strained so that 
both the lattice structure and the local electronic den- 
sity of states (DOS) of graphene are significantly af- 
fected [34]. The non-trivial strain results in signifi- 
cant pseudo-magnetic fields B s up to ~ 50 Tesla, as 
manifested by the DOS peaks at energies correspond- 
ing to quantized Landau levels \E n \ oc -\/n, where 
n denotes both integers and fractional numbers at 0, 
±l,+2,±3,4,6,±l/3,±2/3,±5/3. In addition, signif- 
icant charging effects are found along the strongest 
strained areas of the sample, consistent with strain- 
induced scalar potential. On the other hand, for CVD- 
grown graphene transferred from copper to Si02, both 
the pseudo-magnetic fields and charging effects dimin- 
ish in most parts of the sample except for small re- 
gions containing topological ridges that cannot be re- 
laxed by changing the substrate. The occurrence of 
conductance peaks at quantized integer and fractional 
Landau levels in strained graphene is analogous to the 
magnetic field-induced integer and fractional quantum 
Hall effects (IQHE and FQHE), suggesting that the 
two-dimensional Dirac electrons may be quantum con- 
fined into correlated many-body states by strain. The 
strained-induced charging effects and quantized states 
also enable new possibilities of controlling the energy 
gaps and doping levels of graphene through nano-scale 
"strain engineering". 



2. Experimental 

Our experimental approach to investigating the local 
electronic and structural correlations of graphene on di- 
electric or metallic substrates is to perform STM/STS 
studies using a home built cryogenic STM, which was 
compatible with high magnetic fields and also capable 
of variable temperature control from room temperature 
to 6 K, with a vacuum level of ~ 10~ 10 Torr at the low- 
est temperatures. For studies reported in this work, the 
measurement conditions were at 77 K under high vac- 
uum (< 10 Torr) and in zero magnetic field. The STM 
tips were made of Pt/Ir alloy and were produced by me- 
chanical cleavage and electrochemical polishing. The 
tips thus prepared were first characterized on graphite 
for atomic sharpness and then cleaned by gently touch- 
ing the STM tips to an amorphous gold surface. Both 
topographic and spectroscopic measurements were per- 
formed simultaneously at every location in a (128 x 128) 
pixel grid. At each pixel location, the tunnel junction 
was independently established so that the junction re- 
sistance of 1.5 GQ. was maintained across the sample. 

The differential conductance, {dl/dV), was calcu- 
lated from the best polynomial fit of each tunneling cur- 
rent (/) vs. bias voltage (V) curve. The data acqui- 
sition typically involved averaging ten I-vs.-V spectra 
per pixel. Polynomial fitting to successive segments 
of the averaged I-vs.-V curve at each pixel was made 
over a finite voltage range, and the (dI/dV)-vs.-V spec- 
trum per pixel was obtained by collecting the derivatives 
of the best fitted segments at all voltages. The spec- 
tra thus derived are generally consistent with those ob- 
tained by superposing a modulated bias voltage with a 
lock-in amplifier on the I-vs.-V measurements. How- 
ever, the modulated-bias technique generally requires 
much longer time to collect both the spectroscopic and 
topographical data pixel by pixel over extended areas 
with high spatial resolution, and so was not employed 
in this work. 

The procedures for CVD growth of graphene on 
copper foils were similar to those described in litera- 
ture 13311 . Briefly, graphene films were primarily grown 
on 25 fim thick, highly polished polycrystalline cop- 
per foils in a furnace consisting of a fused silica tube 
heated in a split tube furnace. The fused silica tube 
loaded with copper foils was first evacuated and then 
back filled with pure hydrogen gas, heated to ~ 1000°C, 
and maintained under partial hydrogen pressure with ar- 
gon. A flow of methane was subsequently introduced 
for a desired period of time at a total pressure of 500 
mTorr. Finally, the furnace was cooled to room tem- 
perature, and the copper foils coated with graphene ap- 
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peared brighter relative to the as-received copper foils, 
consistent with previous reports i33ll . The graphene 
sheets thus prepared were found to be largely single- 
or double-layered based on Raman spectroscopic stud- 
ies. Before the STM measurements, a graphene-coated 
copper foil was first cleaned and then loaded onto the 
STM probe. The STM probe was then evacuated to 
high vacuum condition (< 10~ 7 Torr). Two different 
CVD-grown graphene samples on copper were studied 
at 77 K. Both topographic and spectroscopic measure- 
ments on single-layer graphene regions were carried out 
on both samples from relatively large area surveys down 
to atomic-scale investigations. 

To transfer CVD-grown graphene samples to SiC<2 
substrates, SiC<2 substrates were prepared by first ther- 
mally grown a 290 nm thick SiC<2 layer on the sili- 
con wafer (p-type (100)), followed by gentle sonica- 
tion of the substrate in acetone and then pure alcohol 
for about two minutes. The substrate was then baked 
at 1 15°C on a hotplate, nitrogen blown dry while being 
cooled down. For the CVD-grown graphene on cop- 
per, a layer of PMMA was first deposited on top of 
the graphene sample as scaffolding, and then the cop- 
per substrate was removed with nitric acid. Next, the 
PMMA/graphene sample was placed on a SiC<2 sub- 
strate with the graphene side down. Finally, the PMMA 
was removed with acetone. The transferred graphene 
was cleaned and annealed in an argon atmosphere at 
400°C for 30 minutes. Finally, electrodes were created 
by thermally evaporating 2.5 nm chromium and 37.5nm 
gold through an aluminum foil shadow mask. Prior to 
STM measurements, the sample was rinsed again gently 
with acetone followed by pure alcohol to remove possi- 
ble organic surface contaminants. The STM/STS stud- 
ies on one of the transferred graphene samples on SiC<2 
were all conducted on the single-layer region. 



3. Results and Analysis 

Systematic studies of the CVD grown graphene on 
copper revealed large height variations everywhere, as 
exemplified by the topographic images Figs, la - lc and 
the corresponding height histograms in Figs. Id - If for 
one of the two samples. The ripple-like height varia- 
tions are consistent with previous reports 113311 and are 
primarily associated with two physical causes. The first 
is the large difference in the thermal contraction coef- 
ficients of graphene and copper. Namely, upon cooling 
the sample from ~ 1000°C to room temperature after 
the CVD growth of graphene, the expanded graphene 
(due to its negative thermal contraction coefficient) can 
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Figure 1: (Color online) Topographic images of a CVD-grown 
graphene-on-copper sample for (a) a (2300 X 2300) nm 2 area, (b) a 
(50 X 50) nm 2 area and (c) a (3.0 X 3.0) nm 2 area. The white box in 
(a) represents the region that is enlarged into the image in (b), and the 
white box in (b) represents the area that is enlarged into the image in 
(c). The corresponding height histograms for the images shown in (a) 
- (c) are given in (d) - (f). 



no longer remain fiat on the thermally contracted cop- 
per foil, thus resulting in large ripples. The second cause 
is the original surface roughness of the copper foils, as 
manifested in Fig. 2, where the surface morphology im- 
aged by an atomic force microscope (AFM) and the cor- 
responding height histograms of a copper foil processed 
under the same conditions as those for the CVD-grown 
graphene on copper samples are shown is Figs. 2a - 
2d. Comparing Figs. 1 and 2 suggesting that the long- 
wavelength height variations (Az) of the CVD-grown 
graphene on copper may be associated with the surface 
morphology of the copper foil, and the averaged varia- 
tion over a length scale L is estimated at (Az/L) < 0.01. 
On the other hand, significant short-wavelength height 
variations, such as those exemplified in Figs. 2e - 2f over 
a small (3.0 x 3.0) nm 2 area may be primarily attributed 
to the difference in the thermal contraction coefficients, 
where the height variations could reach as much as 1 
nm over a (3.0 x 3.0) nm 2 area for both samples so that 
(Az/L) ~ 0.3, which are much larger than those asso- 
ciated with the long-wavelength height variations and 
therefore suggest significant strain on graphene. 

The large strain associated with the CVD-grown 
graphene-on-copper samples apparently incurred signif- 
icant lattice distortion to the arrangements of carbon 
atoms. As exemplified in the atomically resolved topo- 
graphic image in Fig. 3a, different atomic arrangements 
were apparent in different regions of the (3.0 x 3.0) nm 2 
area in view. Specifically, we found that in a more re- 
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laxed region of the graphene sample (denoted as the a- 
region) the atomic structure appeared to resemble the 
honeycomb or hexagon-like lattice with slight distor- 
tion. On the other hand, for spectra taken on the ar- 
eas slightly below the ridge (j6-region) of a ripple-like 
feature, the atomic arrangements deviated strongly from 
those of graphene, showing either nearly square-lattice 
or disordered atomic structure in the /3-region. Addi- 
tionally, we found certain areas that were completely 
disordered without atomic resolution, as exemplified in 
Fig. 3c and denoted as the y-region. The disordered re- 
gion might be associated with amorphous carbon grown 
on defected features (e.g. the bright features shown in 
the AFM images in Figs. 2a and 2b) on the copper sub- 
strate. 

To investigate the effect of lattice distortion on the lo- 
cal electronic DOS of graphene, we compared the typ- 
ical tunneling conductance (dl/dV) vs. V spectra ob- 
tained from different regions of atomic arrangements. 
As shown in Figs. 3b and 3d, the representative spectra 
for the more relaxed ff-region of the graphene sample 
shown in Fig. 1 resembled those found in mechanically 
exfoliated graphene except a substantially larger zero- 
bias conductance and additional conductance peaks. On 
the other hand, the typical spectra associated with the 
strongly strained /3-region appeared to be asymmetric 
relative to those in the a region, with sharp conductance 
peaks more closely spaced. Additionally, for disordered 
y-regions without atomic resolution, nearly parabolic 
spectra that deviate fundamentally from the linear en- 
ergy dependence of the Dirac fermions were found, as 
exemplified in Fig. 3d. For a- and /3-type spectra exem- 
plified in Fig. 3b, the energy interval between consecu- 
tive conductance peaks appear to decrease with increas- 
ing energy. On the other hand, for y-type spectra, the 
conductance peaks seem randomly spaced. Similar re- 
sults are confirmed in the second CVD-grown graphene- 
on-copper sample, as exemplified in Figs. 3e - 3f, where 
the conductance peaks of the representative point spec- 
tra correspond to a pseudo-magnetic field B s = 35 + 5 
Tesla. 

To better understand the significance of the conduc- 
tance peaks, we subtract off a parabolic conductance 
background associated with the copper contribution (in- 
dicated by the dashed lines in Figs. 3b and 3f) and plot 
the resulting conductance against (E - £birac) where 
£birac is the Dirac energy. As exemplified in Figs. 4a 
and 4c for the point spectra after background subtrac- 
tion, we find distinct peaks occurring at energies pro- 
portional to VN for n — 0, ± 1 , +2, ±3,4 and 6. Further- 
more, for stronger strained /8-type spectra, sharp peaks 
at energies proportional to fractional numbers + y/T/3, 




Figure 2: (Color online) AFM studies of the surface morphology of 
a copper foil processed in the same way as the CVD-grown graphene 
on copper is compared with the STM surface topography of a CVD- 
grown graphene on copper: (a) An AFM image of the surface mor- 
phology of the copper foil over a (5.0x5.0)/jm 2 area, (b) An AFM im- 
age of the surface morphology of the copper foil over a (500x500)nm _ 
area. (c) The height histogram of the surface morphology along the 
red line in (a), (d) The height histogram of the surface morphology 
along the red line in (b). From (a) - (d) we find relatively long- 
wavelength height variations on the order of (Az/L) < 0.01, where 
L denotes the lateral length for total height variation Az. (e) The 
height variations along the dashed line in the STM topography in 
Fig. lc for a CVD-grown graphene on copper, showing atomic steps 
(solid lines) superposed on top of a gradual background (dotted lines), 
(f) The height variations in (e) after removing the long-wavelength 
background, showing sharp steps of either ~ 0.3 nm or ~ 0.6 nm in 
height, comparable to one or two c-axis lattice constants of graphite. 
From (e)-(f) we find short-wavelength height variations on the order of 
(Az/L) ~ 0.3, much larger than the long wavelength variations shown 
in (a) -(d). 
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Figure 3: (Color online) Lattice distortion of the CVD-grown 
graphene-on-copper samples and the resulting effect on the local elec- 
tronic density of states: (a) A representative atomically resolved to- 
pography of one of the samples over a (3.0 X 3.0) nm 2 area, showing 
different atomic arrangements of the graphene sample in view, with 
slightly skewed honeycomb (or hexagonal) lattices in the ff-region and 
a combination of nearly squared and distorted atomic arrangements in 
the /3-region. (b) Representative tunneling conductance dl /dV-vs.-V 
point spectra in the a-region (point 1 ) and /3-region (points 1 and 2) 
shown in (a). The spectra are slighted shifted relative to each other 
along the dl/dV axis for clarify, (c) A representative topography of 
the same sample in (a) over a disordered (3.0 X 3.0) nm 2 area, (d) 
Representative dl/dV-vs.-V point spectra for the specific points in the 
y-region shown in (c). (e) An atomically resolved topography of an- 
other CVD-grown graphene-on-copper sample over a (5.0 X 5.0) nm 2 
area, showing distorted lattices, (f) Representative dl/dV-vs.-V point 
spectra of the sample shown in (e). The spectra for points 1 and 2 
are shifted up slightly along the dl/dV axis for clarity. All data were 
taken at 7/ = 77 K. 
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Figure 4: (Color online) Manifestation of strained-induced pseudo- 
magnetic fields and quantized Landau levels in two CVD-grown 
graphene-on-copper samples at T = 77 K: (a) The tunneling con- 
ductance (dl/dV) spectra vs. (E - fbirac) after the subtraction of 
the parabolic conductance background shown in Fig. 3b. Here £birac 
denotes the Dirac energy. Distinct peaks at energies proportional to 
± VM are clearly visible for n = 0, ±1, ±2, ±3. Additional humps at 
fractional values n = -5/3, -2/3 and ±1/3 are also identifiable, (b) 
(E-E£,i rdc .)-vs- VM for the representative point spectra in (a), showing 
that the data largely fall on linear curves that correspond to pseudo- 
magnetic fields ranging from B s = 29 ± 8 Tesla for the less strained 
a-region to 50 ± 5 Tesla for the stronger strained /3-region. (c) The 
(dl/dV) spectra vs. (E - /ioirac) after the subtraction of the parabolic 
conductance background shown in Fig. 3f. (d) (E - £Dirac) _vs - VM 
for the representative point spectra in (c), showing a pseudo-magnetic 
field B s = 35 ± 5 Tesla. 

- V2/3 and V5/3 are clearly visible, in contrast to the 
weak "humps" occurring in the weaker strained a-type 
spectra at energies proportional to VI /3 and - V5/3. 

By plotting (E - £birac) versus VN data for multiple 
spectra taken in Figs. 4a and 4c, we find the (E - Z?Dirac) 
data taken for each point spectrum in the atomically re- 
solved strained regions all fall approximately on one lin- 
ear curve, as exemplified in Figs. 4b and 4d, and the 
slope for each (E - £'Dirac)-vs-|n| curve increases with 
increasing strain. 

The results manifested in Fig. 4 can be understood 
in terms of the presence of strain-induced pseudo- 
magnetic fields B s , with larger strain giving rise to larger 
B s values ITll[l2lElll . Specifically, the pseudo-Landau 
levels E n of Dirac electrons under a given B s satisfy the 
following relation: 

E n = sga(n)J(2e^ F hB s )\n\, (1) 



where Vf denotes the Fermi velocity of graphene, B s — 
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Figure 6: (Color online) Spatial maps of the effective compres- 
sion/dilation strain components for a (3 X 3) nm 2 area of the CVD- 
grown graphene-on-copper sample shown in Figure 5, where the strain 
components are estimated by considering the spatial variations in the 
height z(x, y): (a) (dz/dx) 2 , (b) (dz/dyf, and (c) [(dz/dx) 2 +(dz/dy) 2 ]. 



Figure 5: (Color online) Constant- voltage tunneling conductance 
maps at quantized energies of the pseudo-Landau levels over the same 
(3 X 3) nm 2 area and for T = 11 K: (a) n = 0, (b) n = 1, (c) 
n = 2, (d) n = 3, (e) n = 1/3, (f) n = 2/3, where the n values 
are referenced to the strain-induced Landau levels for the /Jl point 
spectrum. For smaller n values, an approximately one-dimensional 
high-conductance "line" appears near the topographical ridge where 
the most abrupt changes in height occur, suggesting significant charg- 
ing effects. The confinement of high-conductance region diminishes 
with increasing n. (g) Spatial map of the Dirac-point tunneling con- 
ductance for (dI/dV)(x,y) taken at spatially varying Dirac energy 
V = (E Dk ,Je). 



|VxA|, and the strain-induced gauge field A = A x x+A y y 
is given by 112 ] 

.It 



A x 



(2) 



In Eq. © p = -(Slnf/dlna) ~ 2, where t * 3 eV 
denotes the nearest-neighbor hopping constant, and a 
is the nearest carbon-carbon separation. Additionally, 
the tensor components uy = x,y) are defined as 
Ujj = (duj/di), where Uj is the y'th component of the 
displacement field u = u x x + u y y. 

Using \'f ~ 10 6 m/sec, we find that the pseudo- 
magnetic fields are approximately 29 ± 8 Tesla for the 
a region and 50 + 5 Tesla for the /? region. These 
pseudo-magnetic fields are significant even though 
they are smaller than those found in graphene nano- 
bubbles 12111 . probably due to the spatial inhomogeneity 
of the pseudo-magnetic fields that leads to partial can- 
cellations. Further, the magnetic lengths £b associated 
with the pseudo-magnetic fields are estimated to range 
from 3.5 to 5.5 nm according to the following formula: 



p-2 



(3) 



where cDo is the flux quantum. 

The strain found in the CVD-grown graphene on 
copper is not purely shear but also contains compres- 
sion/dilation components. The latter is theoretically 
predicted to gives rise to an effective scalar potential 



V(x, y) and therefore a static charging effect 111 2l l36l 
3711 in addition to the aforementioned pseudo-magnetic 
field. Here V(x, y) is given by E [H, H : 



V(x,y) = Vo(d x u x + dyUy) = V u, 



(4) 



where Vo ~ 3 eV, and u is the dilation/compression 
strain. For strain induced by corrugations along the 
out of plane direction, the typical strain is of order 
u ~ (Az/L) 2 , where Az denotes the height fluctuations 
and L is the size of the strained region. Although the 
charging effect associated with the scalar potential can 
be suppressed in a graphene layer suspended over a 
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metal if Az <k €b II 1 1 1 Il2ll . the mean value of Az over 
the (3 x 3) nm 2 area of the CVD-grown graphene on 
copper was Az ~ 1 nm, comparable to Cg. This quan- 
titative consideration therefore suggests that significant 
charging effects may be observed. 

To visualize the charging effect, we show in Figs. 5a 
- 5f the conductance maps for one of the CVD- 
grown graphene-on-copper samples at different con- 
stant bias voltages over the same (3 x 3) nm 2 area. 
The bias voltages are chosen so that they correspond 
to the pseudo-Landau levels of the f3\ point spectrum 
for n = 0, 1, 2, 3, 1/3 and 2/3. A nearly one- 
dimensional high-conductance region following closely 
to the most strained region immediately below the ridge 
is clearly visible in maps associated with smaller n val- 
ues, confirming the notion of strain-induced charging 
effects. For higher energies (larger n values), the high- 
conductance region becomes less confined, which is 
reasonable because of the higher confinement energies 
(oc -\Jn) required for the Dirac electrons. We further 
note that the strain-induced charging effect also results 
in shifts in the Dirac energy. As manifested by the 
spatial map in Fig. 5g for the tunneling conductance 
taken at the spatially varying Dirac energy Euimc(.x,y) 
for the same graphene sample, the {dI/dV)(x,y) map 
for V = (£birac/e) also exhibits excess conductance in 
strongly strained regions. 
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Figure 7: (Color online) Topographic images of a CVD-grown 
graphene transferred from copper to SiC>2 for (a) a (100 X 100) nm 2 
area and (b) a (10 X 10) nm 2 area. Both images were taken at 77 K. 
The corresponding height histograms for the images shown in (a) and 
(b) are given in (c) and (d), respectively. These height variations are 
much smaller than those found in Fig. la - lc. 



To further manifest the compression/dilation strain- 
induced charging effect, we show in Figs. 6a - 6c spa- 
tial maps of effective compression/dilation strain com- 
ponents (dz/dx) 2 , (dz/dy) 2 and [(dz/dx) 2 + (dz/dy) 2 ]. 
Comparing Fig. 6c with the conductance maps in Fig. 5, 
particularly Fig. 5g, we find apparent correlation be- 
tween the magnitude of strain and the conductance at 
the Dirac energy. 

Next, we consider comparative studies of a CVD- 
grown graphene sample transferred from copper to 
SiC>2. As shown in Figs. 7a - 7d, the height variations 
in the transferred graphene sample are much reduced 
relative to those shown in Figs, la - lc. The signif- 
icant reduction in the average height variations of the 
ripples results in lattice structures and tunneling con- 
ductance spectra of the graphene sample similar to those 
found in mechanically exfoliated graphene on SiC»2 |@], 
as exemplified in Figs. 8a and 8b. Moreover, the tunnel- 
ing conductance spectra appear to be much smoother 
with no discernible conductance peaks in most of the 
transferred graphene-on-SiC<2 sample, as exemplified by 
the spectra in Figs. 8b and 8e for the points marked in 
Fig. 8a. On the other hand, occasional tunneling spec- 
tra near strained regions of the transferred graphene still 
revealed conductance peaks at quantized energies, as 
exemplified by the point spectrum in Figs. 8c and 8f. 
These regions contain topological ridges that cannot be 
relaxed by transferring the graphene sample from cop- 



per to SiC<2 substrate, so that the corresponding pseudo- 
magnetic field remains substantial, as exemplified in 
Fig. 8d where B s =16 + 2 Tesla. 

Similarly, the reduced strain for the transferred 
graphene-on-SiC>2 sample leads to reduced charging 
effect, as manifested by the overall lower and more 
spatially homogeneous Dirac -point conductance of the 
transferred graphene sample in Fig. 9c relative to that of 
the CVD-grown graphene-on-copper sample in Fig. 9a. 
Further comparison of the effective lattice distortion 
[(dz/dx) 2 + (dz/dy) 2 ] between the transferred graphene 
and the graphene-on-copper in Figs. 9b and 9d reveal 
that the compression/dilation strain is much reduced in 
the transferred sample. For strongly strained graphene, 
the magnitude of the compression/dilation strain corre- 
lates well with the Dirac -point conductance, as mani- 
fested by Figs. 9a and 9b. On the other hand, for negli- 
gible compression/dilation strain as in the case of trans- 
ferred graphene-on-Si02, there is no longer direct corre- 
lation between the Dirac -point conductance and the lo- 
cal strain, because the former may be largely dependent 
on the local charge impurities on substrates. Overall, 
our STM/STS studies of various graphene samples have 
manifested the strong influence of strain and substrates 
on the structural rippling effect and the electronic DOS 
of graphene. 

4. Discussion 

We have shown that significant differences between 
the thermal contraction coefficients of graphene and the 
transition-metal substrates for CVD growth graphene 
result in strong and non-uniform lattice distortion at the 
nano-scale, which provides a new ground for investi- 
gating the strain-induced charging and pseudo-magnetic 
field effects. In particular, the strain-induced giant 
pseudo-magnetic fields enable studies of the integer 
and fractional quantum Hall effects associated with the 
pseudo-spin (i.e., the valley) degrees of freedom by ex- 
ploiting the unique strength of the pseudo-spin-orbit 
coupling in graphene. This realization is analogous to 
the situation of two-dimensional topological insulators 
where a gapped bulk state coexists with the counter- 
propagating states at the boundaries without breaking 
the time-reversal symmetry I13lll4lll5ll . While in topo- 
logical insulators the spin-orbit interaction is generally 
sizable (10~' ~ 1 eV), in graphene the spin-orbit cou- 
pling is very small so that the corresponding Landau 
gaps are smaller than 1 ji/eV. Therefore, direct manifes- 
tation of the quantum spin Hall effect in graphene can 
only be realized under extreme empirical conditions of 
ultra-low temperatures (< 10 mK) and ultra-high carrier 
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Figure 8: (Color online) Topography and spectroscopy taken at 77 K on CVD-grown graphene transferred to Si02: (a) A representative topography 
over a (10 X 10) nm 2 area, showing gradual height variations and therefore corresponding to a relaxed region, (b) Representative tunneling 
conductance spectra of the relaxed region at multiple points marked in (a). The zero-bias conductance for the spectra taken at points a and b has 
been offset for easier comparison. Clearly all spectra are smooth without conductance peaks associated with pseudo-magnetic fields, in contrast 
to the findings for CVD-grown graphene-on-copper samples shown in Figs. 3 and 4. (c) A tunneling conductance spectrum in the transferred 
graphene-on-Si02 sample that corresponds a strained region (topography not shown). The spectrum reveals distinct quantized conductance peaks 
spaced at much smaller energies than those observed in Figure 4 for the CVD-grown graphene-on-copper samples, (d) (E - £Dirac) _vs -- VM ror 
the spectrum taken in the strained region, showing that the conductance peaks fall on one linear curve that corresponds to B s = 16 ± 2 Tesla. (e) 
Multiple tunneling conductance spectra of the relaxed region subtracted by a background, showing no discernible conductance peaks and therefore 
implying B s — > 0. (f) Tunneling conductance spectrum of the strained region subtracted by a background (indicated by the dashed line in (c)), 
showing conductance peaks at quantized energies, with \n\ = 1/3, 1, 2 and 3. 
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Figure 9: (Color online) Strain effects on the Dirac-point tunneling 
conductance at 7/ = 77 K: (a) A (3 X 3) nm 2 spatial map of the Dirac 
point (E = EDirac) tunneling conductance of the same CVD-grown 
graphene-on-copper sample as in Figure 3a. (b) Spatial map of the 
compression/dilation strain [(dzldx) 2 + (dz/dy) 2 ] over the same sam- 
ple area as in (a). The maps in (a) and (b) are clearly well corre- 
lated, (c) A spatial map of the tunneling conductance of the trans- 
ferred graphene-on-Si02 sample taken at the Dirac point i?Dirac over 
a (3 X 3) nm" area, (d) The spatial map of the effective compres- 
sion/dilation strain [(dz/dx) 2 + (dz/dy) 2 ] over the same sample area as 
in (c). The effective strain in (d) is several orders of magnitude smaller 
than that in (b). Hence, the lack of correlation between the maps in 
(c) and (d) implies that the charging effect in relaxed graphene is not 
directly due to strain and may be attributed to charge impurities on the 
substrate. 

mobilities (> 10 7 cm 2 V s _1 ). In contrast, the pseudo- 
magnetic field-induced quantum Hall effect associated 
with the pseudo-spins is easily identifiable in strained 
graphene, as demonstrated by this work and by similar 
studies of graphene nano- bubbles on Pt(lll) DM]. 

It is worth pointing out that our observation of 
pseudo-magnetic field-induced fractionally filled Lan- 
dau levels in graphene may be attributed to the compli- 
cated strain distribution at the nano-scale, which ensures 
both large pseudo-magnetic fields and ballistic behav- 
ior of the Dirac electrons. Hence, strain-induced charge 
confinement in two dimensions (similar to the exter- 
nal magnetic field-induced charge confinement) gives 
rise to correlated many-body interactions. Moreover, 
the significantly more complicated strain distribution 
found in the CVD-grown graphene-on-copper differs 
from the controlled strain in the case of graphene nano- 
bubbles ll2lll where the effective gauge potential leads 
to a relatively uniform magnetic field over each nano- 
bubble. Therefore, further corrections to the Hamilto- 



nian of the Dirac electrons are necessary in the case 
of CVD-grown graphene-on-copper, which may lead to 
additional terms in the effective gauge potential, simi- 
lar to the presence of the Chern-Simons gauge potential 
terms in the case of fractional quantum Hall (FQH) phe- 
nomena. However, self-consistent theoretical studies 
will be necessary to account for the observation of FQH- 
like effects in the CVD-grown graphene-on-copper. 

We further remark that the tunneling spectra acquired 
from highly strained CVD-grown graphene-on-copper 
samples contain not only spectral characteristics of the 
strained graphene but also excess conduction channels 
associated with the underlying copper substrate. The 
excess conduction channels may be responsible for the 
commonly observed suppression of the zero-bias DOS 
associated with the n — Landau level. However, the 
general effects of the underlying copper substrate on the 
graphene spectra are still unknown and await further 
theoretical investigation. 

The findings reported in this work are highly rel- 
evant to the development of graphene-based nano- 
electronics. In particular, proper "strain engineering" 
may be achieved by designing the substrate for graphene 
so that either local doping or energy-gap engineering as- 
sociated with the pseudo-Landau levels can be realized 
at the nano-scale. Ultimately, strain-induced charging 
effects and pseudo-magnetic fields may also be achieved 
at the macroscopic scale if quasi-periodic occurrences 
of these nano-scale effects can be constructed. 

5. Conclusions 

We have conducted spatially resolved topographic 
and spectroscopic studies of CVD-grown graphene-on- 
copper and transferred graphene-on-SiCh samples. Our 
investigation reveals the important influence of the sub- 
strate and strain on the carbon atomic arrangements 
and the electronic DOS of graphene. For CVD-grown 
graphene remaining on the copper substrate, the mono- 
layer carbon structures exhibit ripples and are strongly 
strained, with different regions exhibiting varying lat- 
tice structures and electronic density of states. In par- 
ticular, topographical ridges appear along the bound- 
aries between different lattice structures, which also ex- 
hibit excess charging effects. Additionally, the large 
and non-uniform strain induces pseudo-magnetic fields 
up to ~ 50 Tesla, as manifested by the tunneling con- 
ductance peaks at quantized energies that are associated 
with both integer and fractional pseudo-magnetic field- 
induced Landau levels. In contrast, for graphene trans- 
ferred from copper to S1O2 after CVD growth, the typ- 
ical graphene structure is largely restored and the aver- 



9 



age strain on the whole is much reduced, so are the cor- 
responding charging effects and pseudo-magnetic fields, 
with the exception of limited areas that are associ- 
ated with topological defects and therefore the result- 
ing strain cannot be relaxed after the transfer of sample 
from copper to S1O2 substrates. Our findings suggest 
feasible nano-scale strain engineering of the electronic 
states of CVD-grown graphene by proper design of the 
substrates and growth conditions. 
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